3.462 \(\int \frac {(c+d \sin (e+f x))^4}{(a+a \sin (e+f x))^2} \, dx\)

Optimal. Leaf size=195 \[ \frac {d^2 \left (2 c^2+16 c d-21 d^2\right ) \sin (e+f x) \cos (e+f x)}{6 a^2 f}+\frac {d^2 x \left (12 c^2-16 c d+7 d^2\right )}{2 a^2}+\frac {2 d \left (c^3+8 c^2 d-20 c d^2+8 d^3\right ) \cos (e+f x)}{3 a^2 f}-\frac {(c-d) (c+8 d) \cos (e+f x) (c+d \sin (e+f x))^2}{3 a^2 f (\sin (e+f x)+1)}-\frac {(c-d) \cos (e+f x) (c+d \sin (e+f x))^3}{3 f (a \sin (e+f x)+a)^2} \]

[Out]

1/2*d^2*(12*c^2-16*c*d+7*d^2)*x/a^2+2/3*d*(c^3+8*c^2*d-20*c*d^2+8*d^3)*cos(f*x+e)/a^2/f+1/6*d^2*(2*c^2+16*c*d-
21*d^2)*cos(f*x+e)*sin(f*x+e)/a^2/f-1/3*(c-d)*(c+8*d)*cos(f*x+e)*(c+d*sin(f*x+e))^2/a^2/f/(1+sin(f*x+e))-1/3*(
c-d)*cos(f*x+e)*(c+d*sin(f*x+e))^3/f/(a+a*sin(f*x+e))^2

________________________________________________________________________________________

Rubi [A]  time = 0.36, antiderivative size = 195, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, integrand size = 25, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.120, Rules used = {2765, 2977, 2734} \[ \frac {2 d \left (8 c^2 d+c^3-20 c d^2+8 d^3\right ) \cos (e+f x)}{3 a^2 f}+\frac {d^2 \left (2 c^2+16 c d-21 d^2\right ) \sin (e+f x) \cos (e+f x)}{6 a^2 f}+\frac {d^2 x \left (12 c^2-16 c d+7 d^2\right )}{2 a^2}-\frac {(c-d) (c+8 d) \cos (e+f x) (c+d \sin (e+f x))^2}{3 a^2 f (\sin (e+f x)+1)}-\frac {(c-d) \cos (e+f x) (c+d \sin (e+f x))^3}{3 f (a \sin (e+f x)+a)^2} \]

Antiderivative was successfully verified.

[In]

Int[(c + d*Sin[e + f*x])^4/(a + a*Sin[e + f*x])^2,x]

[Out]

(d^2*(12*c^2 - 16*c*d + 7*d^2)*x)/(2*a^2) + (2*d*(c^3 + 8*c^2*d - 20*c*d^2 + 8*d^3)*Cos[e + f*x])/(3*a^2*f) +
(d^2*(2*c^2 + 16*c*d - 21*d^2)*Cos[e + f*x]*Sin[e + f*x])/(6*a^2*f) - ((c - d)*(c + 8*d)*Cos[e + f*x]*(c + d*S
in[e + f*x])^2)/(3*a^2*f*(1 + Sin[e + f*x])) - ((c - d)*Cos[e + f*x]*(c + d*Sin[e + f*x])^3)/(3*f*(a + a*Sin[e
 + f*x])^2)

Rule 2734

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Simp[((2*a*c
+ b*d)*x)/2, x] + (-Simp[((b*c + a*d)*Cos[e + f*x])/f, x] - Simp[(b*d*Cos[e + f*x]*Sin[e + f*x])/(2*f), x]) /;
 FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0]

Rule 2765

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Sim
p[((b*c - a*d)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*(c + d*Sin[e + f*x])^(n - 1))/(a*f*(2*m + 1)), x] + Dist[1/
(a*b*(2*m + 1)), Int[(a + b*Sin[e + f*x])^(m + 1)*(c + d*Sin[e + f*x])^(n - 2)*Simp[b*(c^2*(m + 1) + d^2*(n -
1)) + a*c*d*(m - n + 1) + d*(a*d*(m - n + 1) + b*c*(m + n))*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, d, e,
f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && LtQ[m, -1] && GtQ[n, 1] && (IntegersQ
[2*m, 2*n] || (IntegerQ[m] && EqQ[c, 0]))

Rule 2977

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[((A*b - a*B)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*(c + d*Sin[e + f*x]
)^n)/(a*f*(2*m + 1)), x] - Dist[1/(a*b*(2*m + 1)), Int[(a + b*Sin[e + f*x])^(m + 1)*(c + d*Sin[e + f*x])^(n -
1)*Simp[A*(a*d*n - b*c*(m + 1)) - B*(a*c*m + b*d*n) - d*(a*B*(m - n) + A*b*(m + n + 1))*Sin[e + f*x], x], x],
x] /; FreeQ[{a, b, c, d, e, f, A, B}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && LtQ
[m, -2^(-1)] && GtQ[n, 0] && IntegerQ[2*m] && (IntegerQ[2*n] || EqQ[c, 0])

Rubi steps

\begin {align*} \int \frac {(c+d \sin (e+f x))^4}{(a+a \sin (e+f x))^2} \, dx &=-\frac {(c-d) \cos (e+f x) (c+d \sin (e+f x))^3}{3 f (a+a \sin (e+f x))^2}-\frac {\int \frac {(c+d \sin (e+f x))^2 \left (-a \left (c^2+5 c d-3 d^2\right )+a (2 c-5 d) d \sin (e+f x)\right )}{a+a \sin (e+f x)} \, dx}{3 a^2}\\ &=-\frac {(c-d) (c+8 d) \cos (e+f x) (c+d \sin (e+f x))^2}{3 a^2 f (1+\sin (e+f x))}-\frac {(c-d) \cos (e+f x) (c+d \sin (e+f x))^3}{3 f (a+a \sin (e+f x))^2}-\frac {\int (c+d \sin (e+f x)) \left (-a^2 (19 c-16 d) d^2+a^2 d \left (2 c^2+16 c d-21 d^2\right ) \sin (e+f x)\right ) \, dx}{3 a^4}\\ &=\frac {d^2 \left (12 c^2-16 c d+7 d^2\right ) x}{2 a^2}+\frac {2 d \left (c^3+8 c^2 d-20 c d^2+8 d^3\right ) \cos (e+f x)}{3 a^2 f}+\frac {d^2 \left (2 c^2+16 c d-21 d^2\right ) \cos (e+f x) \sin (e+f x)}{6 a^2 f}-\frac {(c-d) (c+8 d) \cos (e+f x) (c+d \sin (e+f x))^2}{3 a^2 f (1+\sin (e+f x))}-\frac {(c-d) \cos (e+f x) (c+d \sin (e+f x))^3}{3 f (a+a \sin (e+f x))^2}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]  time = 1.95, size = 378, normalized size = 1.94 \[ \frac {\left (\sin \left (\frac {1}{2} (e+f x)\right )+\cos \left (\frac {1}{2} (e+f x)\right )\right ) \left (3 d \cos \left (\frac {1}{2} (e+f x)\right ) \left (64 c^3+48 c^2 d (3 e+3 f x-4)-32 c d^2 (6 e+6 f x-5)+7 d^3 (12 e+12 f x-7)\right )-\cos \left (\frac {3}{2} (e+f x)\right ) \left (16 c^4+128 c^3 d+48 c^2 d^2 (3 e+3 f x-10)-16 c d^3 (12 e+12 f x-41)+d^4 (84 e+84 f x-239)\right )+3 \left (2 \sin \left (\frac {1}{2} (e+f x)\right ) \left (8 c^4+32 c^3 d+d^2 \cos (e+f x) \left (48 c^2 (e+f x)-64 c d (e+f x+1)+d^2 (28 e+28 f x+27)\right )+96 c^2 d^2 e+96 c^2 d^2 f x-144 c^2 d^2-2 d^3 (8 c-3 d) \cos (2 (e+f x))-128 c d^3 e-128 c d^3 f x+144 c d^3+d^4 \cos (3 (e+f x))+56 d^4 e+56 d^4 f x-50 d^4\right )+d^3 (16 c-5 d) \cos \left (\frac {5}{2} (e+f x)\right )+d^4 \cos \left (\frac {7}{2} (e+f x)\right )\right )\right )}{48 a^2 f (\sin (e+f x)+1)^2} \]

Antiderivative was successfully verified.

[In]

Integrate[(c + d*Sin[e + f*x])^4/(a + a*Sin[e + f*x])^2,x]

[Out]

((Cos[(e + f*x)/2] + Sin[(e + f*x)/2])*(3*d*(64*c^3 + 48*c^2*d*(-4 + 3*e + 3*f*x) - 32*c*d^2*(-5 + 6*e + 6*f*x
) + 7*d^3*(-7 + 12*e + 12*f*x))*Cos[(e + f*x)/2] - (16*c^4 + 128*c^3*d + 48*c^2*d^2*(-10 + 3*e + 3*f*x) - 16*c
*d^3*(-41 + 12*e + 12*f*x) + d^4*(-239 + 84*e + 84*f*x))*Cos[(3*(e + f*x))/2] + 3*((16*c - 5*d)*d^3*Cos[(5*(e
+ f*x))/2] + d^4*Cos[(7*(e + f*x))/2] + 2*(8*c^4 + 32*c^3*d - 144*c^2*d^2 + 144*c*d^3 - 50*d^4 + 96*c^2*d^2*e
- 128*c*d^3*e + 56*d^4*e + 96*c^2*d^2*f*x - 128*c*d^3*f*x + 56*d^4*f*x + d^2*(48*c^2*(e + f*x) - 64*c*d*(1 + e
 + f*x) + d^2*(27 + 28*e + 28*f*x))*Cos[e + f*x] - 2*(8*c - 3*d)*d^3*Cos[2*(e + f*x)] + d^4*Cos[3*(e + f*x)])*
Sin[(e + f*x)/2])))/(48*a^2*f*(1 + Sin[e + f*x])^2)

________________________________________________________________________________________

fricas [B]  time = 0.45, size = 440, normalized size = 2.26 \[ -\frac {3 \, d^{4} \cos \left (f x + e\right )^{4} - 2 \, c^{4} + 8 \, c^{3} d - 12 \, c^{2} d^{2} + 8 \, c d^{3} - 2 \, d^{4} + 6 \, {\left (4 \, c d^{3} - d^{4}\right )} \cos \left (f x + e\right )^{3} + 6 \, {\left (12 \, c^{2} d^{2} - 16 \, c d^{3} + 7 \, d^{4}\right )} f x - {\left (2 \, c^{4} + 16 \, c^{3} d - 60 \, c^{2} d^{2} + 88 \, c d^{3} - 31 \, d^{4} + 3 \, {\left (12 \, c^{2} d^{2} - 16 \, c d^{3} + 7 \, d^{4}\right )} f x\right )} \cos \left (f x + e\right )^{2} - {\left (4 \, c^{4} + 8 \, c^{3} d - 48 \, c^{2} d^{2} + 104 \, c d^{3} - 38 \, d^{4} - 3 \, {\left (12 \, c^{2} d^{2} - 16 \, c d^{3} + 7 \, d^{4}\right )} f x\right )} \cos \left (f x + e\right ) + {\left (3 \, d^{4} \cos \left (f x + e\right )^{3} + 2 \, c^{4} - 8 \, c^{3} d + 12 \, c^{2} d^{2} - 8 \, c d^{3} + 2 \, d^{4} + 6 \, {\left (12 \, c^{2} d^{2} - 16 \, c d^{3} + 7 \, d^{4}\right )} f x - 3 \, {\left (8 \, c d^{3} - 3 \, d^{4}\right )} \cos \left (f x + e\right )^{2} - {\left (2 \, c^{4} + 16 \, c^{3} d - 60 \, c^{2} d^{2} + 112 \, c d^{3} - 40 \, d^{4} - 3 \, {\left (12 \, c^{2} d^{2} - 16 \, c d^{3} + 7 \, d^{4}\right )} f x\right )} \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )}{6 \, {\left (a^{2} f \cos \left (f x + e\right )^{2} - a^{2} f \cos \left (f x + e\right ) - 2 \, a^{2} f - {\left (a^{2} f \cos \left (f x + e\right ) + 2 \, a^{2} f\right )} \sin \left (f x + e\right )\right )}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c+d*sin(f*x+e))^4/(a+a*sin(f*x+e))^2,x, algorithm="fricas")

[Out]

-1/6*(3*d^4*cos(f*x + e)^4 - 2*c^4 + 8*c^3*d - 12*c^2*d^2 + 8*c*d^3 - 2*d^4 + 6*(4*c*d^3 - d^4)*cos(f*x + e)^3
 + 6*(12*c^2*d^2 - 16*c*d^3 + 7*d^4)*f*x - (2*c^4 + 16*c^3*d - 60*c^2*d^2 + 88*c*d^3 - 31*d^4 + 3*(12*c^2*d^2
- 16*c*d^3 + 7*d^4)*f*x)*cos(f*x + e)^2 - (4*c^4 + 8*c^3*d - 48*c^2*d^2 + 104*c*d^3 - 38*d^4 - 3*(12*c^2*d^2 -
 16*c*d^3 + 7*d^4)*f*x)*cos(f*x + e) + (3*d^4*cos(f*x + e)^3 + 2*c^4 - 8*c^3*d + 12*c^2*d^2 - 8*c*d^3 + 2*d^4
+ 6*(12*c^2*d^2 - 16*c*d^3 + 7*d^4)*f*x - 3*(8*c*d^3 - 3*d^4)*cos(f*x + e)^2 - (2*c^4 + 16*c^3*d - 60*c^2*d^2
+ 112*c*d^3 - 40*d^4 - 3*(12*c^2*d^2 - 16*c*d^3 + 7*d^4)*f*x)*cos(f*x + e))*sin(f*x + e))/(a^2*f*cos(f*x + e)^
2 - a^2*f*cos(f*x + e) - 2*a^2*f - (a^2*f*cos(f*x + e) + 2*a^2*f)*sin(f*x + e))

________________________________________________________________________________________

giac [A]  time = 0.24, size = 338, normalized size = 1.73 \[ \frac {\frac {3 \, {\left (12 \, c^{2} d^{2} - 16 \, c d^{3} + 7 \, d^{4}\right )} {\left (f x + e\right )}}{a^{2}} + \frac {6 \, {\left (d^{4} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{3} - 8 \, c d^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 4 \, d^{4} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} - d^{4} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 8 \, c d^{3} + 4 \, d^{4}\right )}}{{\left (\tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 1\right )}^{2} a^{2}} - \frac {4 \, {\left (3 \, c^{4} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} - 18 \, c^{2} d^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 24 \, c d^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} - 9 \, d^{4} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 3 \, c^{4} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 12 \, c^{3} d \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 54 \, c^{2} d^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 60 \, c d^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 21 \, d^{4} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 2 \, c^{4} + 4 \, c^{3} d - 24 \, c^{2} d^{2} + 28 \, c d^{3} - 10 \, d^{4}\right )}}{a^{2} {\left (\tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 1\right )}^{3}}}{6 \, f} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c+d*sin(f*x+e))^4/(a+a*sin(f*x+e))^2,x, algorithm="giac")

[Out]

1/6*(3*(12*c^2*d^2 - 16*c*d^3 + 7*d^4)*(f*x + e)/a^2 + 6*(d^4*tan(1/2*f*x + 1/2*e)^3 - 8*c*d^3*tan(1/2*f*x + 1
/2*e)^2 + 4*d^4*tan(1/2*f*x + 1/2*e)^2 - d^4*tan(1/2*f*x + 1/2*e) - 8*c*d^3 + 4*d^4)/((tan(1/2*f*x + 1/2*e)^2
+ 1)^2*a^2) - 4*(3*c^4*tan(1/2*f*x + 1/2*e)^2 - 18*c^2*d^2*tan(1/2*f*x + 1/2*e)^2 + 24*c*d^3*tan(1/2*f*x + 1/2
*e)^2 - 9*d^4*tan(1/2*f*x + 1/2*e)^2 + 3*c^4*tan(1/2*f*x + 1/2*e) + 12*c^3*d*tan(1/2*f*x + 1/2*e) - 54*c^2*d^2
*tan(1/2*f*x + 1/2*e) + 60*c*d^3*tan(1/2*f*x + 1/2*e) - 21*d^4*tan(1/2*f*x + 1/2*e) + 2*c^4 + 4*c^3*d - 24*c^2
*d^2 + 28*c*d^3 - 10*d^4)/(a^2*(tan(1/2*f*x + 1/2*e) + 1)^3))/f

________________________________________________________________________________________

maple [B]  time = 0.30, size = 618, normalized size = 3.17 \[ \frac {d^{4} \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a^{2} f \left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{2}}-\frac {8 d^{3} \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right ) c}{a^{2} f \left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{2}}+\frac {4 d^{4} \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a^{2} f \left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{2}}-\frac {d^{4} \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{a^{2} f \left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{2}}-\frac {8 d^{3} c}{a^{2} f \left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{2}}+\frac {4 d^{4}}{a^{2} f \left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{2}}+\frac {12 d^{2} \arctan \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )\right ) c^{2}}{a^{2} f}-\frac {16 d^{3} \arctan \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )\right ) c}{a^{2} f}+\frac {7 d^{4} \arctan \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{a^{2} f}-\frac {2 c^{4}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )}+\frac {12 c^{2} d^{2}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )}-\frac {16 c \,d^{3}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )}+\frac {6 d^{4}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )}+\frac {2 c^{4}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{2}}-\frac {8 c^{3} d}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{2}}+\frac {12 c^{2} d^{2}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{2}}-\frac {8 c \,d^{3}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{2}}+\frac {2 d^{4}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{2}}-\frac {4 c^{4}}{3 a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{3}}+\frac {16 c^{3} d}{3 a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{3}}-\frac {8 c^{2} d^{2}}{a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{3}}+\frac {16 c \,d^{3}}{3 a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{3}}-\frac {4 d^{4}}{3 a^{2} f \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )+1\right )^{3}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((c+d*sin(f*x+e))^4/(a+a*sin(f*x+e))^2,x)

[Out]

1/a^2/f*d^4/(1+tan(1/2*f*x+1/2*e)^2)^2*tan(1/2*f*x+1/2*e)^3-8/a^2/f*d^3/(1+tan(1/2*f*x+1/2*e)^2)^2*tan(1/2*f*x
+1/2*e)^2*c+4/a^2/f*d^4/(1+tan(1/2*f*x+1/2*e)^2)^2*tan(1/2*f*x+1/2*e)^2-1/a^2/f*d^4/(1+tan(1/2*f*x+1/2*e)^2)^2
*tan(1/2*f*x+1/2*e)-8/a^2/f*d^3/(1+tan(1/2*f*x+1/2*e)^2)^2*c+4/a^2/f*d^4/(1+tan(1/2*f*x+1/2*e)^2)^2+12/a^2/f*d
^2*arctan(tan(1/2*f*x+1/2*e))*c^2-16/a^2/f*d^3*arctan(tan(1/2*f*x+1/2*e))*c+7/a^2/f*d^4*arctan(tan(1/2*f*x+1/2
*e))-2*c^4/a^2/f/(tan(1/2*f*x+1/2*e)+1)+12/a^2/f/(tan(1/2*f*x+1/2*e)+1)*c^2*d^2-16/a^2/f/(tan(1/2*f*x+1/2*e)+1
)*c*d^3+6/a^2/f/(tan(1/2*f*x+1/2*e)+1)*d^4+2*c^4/a^2/f/(tan(1/2*f*x+1/2*e)+1)^2-8/a^2/f/(tan(1/2*f*x+1/2*e)+1)
^2*c^3*d+12/a^2/f/(tan(1/2*f*x+1/2*e)+1)^2*c^2*d^2-8/a^2/f/(tan(1/2*f*x+1/2*e)+1)^2*c*d^3+2/a^2/f/(tan(1/2*f*x
+1/2*e)+1)^2*d^4-4/3*c^4/a^2/f/(tan(1/2*f*x+1/2*e)+1)^3+16/3/a^2/f/(tan(1/2*f*x+1/2*e)+1)^3*c^3*d-8/a^2/f/(tan
(1/2*f*x+1/2*e)+1)^3*c^2*d^2+16/3/a^2/f/(tan(1/2*f*x+1/2*e)+1)^3*c*d^3-4/3/a^2/f/(tan(1/2*f*x+1/2*e)+1)^3*d^4

________________________________________________________________________________________

maxima [B]  time = 0.45, size = 908, normalized size = 4.66 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c+d*sin(f*x+e))^4/(a+a*sin(f*x+e))^2,x, algorithm="maxima")

[Out]

1/3*(d^4*((75*sin(f*x + e)/(cos(f*x + e) + 1) + 97*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 126*sin(f*x + e)^3/(c
os(f*x + e) + 1)^3 + 98*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 63*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 21*sin(
f*x + e)^6/(cos(f*x + e) + 1)^6 + 32)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 5*a^2*sin(f*x + e)^2/(cos
(f*x + e) + 1)^2 + 7*a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 7*a^2*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 5*a
^2*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 3*a^2*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 + a^2*sin(f*x + e)^7/(cos(f
*x + e) + 1)^7) + 21*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/a^2) - 16*c*d^3*((12*sin(f*x + e)/(cos(f*x + e) +
 1) + 11*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 9*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 3*sin(f*x + e)^4/(cos(f
*x + e) + 1)^4 + 5)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 4*a^2*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 +
 4*a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 3*a^2*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + a^2*sin(f*x + e)^5/(c
os(f*x + e) + 1)^5) + 3*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/a^2) + 12*c^2*d^2*((9*sin(f*x + e)/(cos(f*x +
e) + 1) + 3*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 4)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 3*a^2*sin(
f*x + e)^2/(cos(f*x + e) + 1)^2 + a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3) + 3*arctan(sin(f*x + e)/(cos(f*x +
e) + 1))/a^2) - 2*c^4*(3*sin(f*x + e)/(cos(f*x + e) + 1) + 3*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 2)/(a^2 + 3
*a^2*sin(f*x + e)/(cos(f*x + e) + 1) + 3*a^2*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + a^2*sin(f*x + e)^3/(cos(f*x
 + e) + 1)^3) - 8*c^3*d*(3*sin(f*x + e)/(cos(f*x + e) + 1) + 1)/(a^2 + 3*a^2*sin(f*x + e)/(cos(f*x + e) + 1) +
 3*a^2*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + a^2*sin(f*x + e)^3/(cos(f*x + e) + 1)^3))/f

________________________________________________________________________________________

mupad [B]  time = 9.21, size = 478, normalized size = 2.45 \[ \frac {d^2\,\mathrm {atan}\left (\frac {d^2\,\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )\,\left (12\,c^2-16\,c\,d+7\,d^2\right )}{12\,c^2\,d^2-16\,c\,d^3+7\,d^4}\right )\,\left (12\,c^2-16\,c\,d+7\,d^2\right )}{a^2\,f}-\frac {{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^5\,\left (2\,c^4+8\,c^3\,d-36\,c^2\,d^2+48\,c\,d^3-21\,d^4\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^3\,\left (4\,c^4+16\,c^3\,d-72\,c^2\,d^2+112\,c\,d^3-42\,d^4\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^4\,\left (\frac {16\,c^4}{3}+\frac {8\,c^3\,d}{3}-40\,c^2\,d^2+\frac {224\,c\,d^3}{3}-\frac {98\,d^4}{3}\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^2\,\left (\frac {14\,c^4}{3}+\frac {16\,c^3\,d}{3}-44\,c^2\,d^2+\frac {256\,c\,d^3}{3}-\frac {97\,d^4}{3}\right )+\frac {80\,c\,d^3}{3}+\frac {8\,c^3\,d}{3}+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^6\,\left (2\,c^4-12\,c^2\,d^2+16\,c\,d^3-7\,d^4\right )+\frac {4\,c^4}{3}-\frac {32\,d^4}{3}+\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )\,\left (2\,c^4+8\,c^3\,d-36\,c^2\,d^2+64\,c\,d^3-25\,d^4\right )-16\,c^2\,d^2}{f\,\left (a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^7+3\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^6+5\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^5+7\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^4+7\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^3+5\,a^2\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^2+3\,a^2\,\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )+a^2\right )} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((c + d*sin(e + f*x))^4/(a + a*sin(e + f*x))^2,x)

[Out]

(d^2*atan((d^2*tan(e/2 + (f*x)/2)*(12*c^2 - 16*c*d + 7*d^2))/(7*d^4 - 16*c*d^3 + 12*c^2*d^2))*(12*c^2 - 16*c*d
 + 7*d^2))/(a^2*f) - (tan(e/2 + (f*x)/2)^5*(48*c*d^3 + 8*c^3*d + 2*c^4 - 21*d^4 - 36*c^2*d^2) + tan(e/2 + (f*x
)/2)^3*(112*c*d^3 + 16*c^3*d + 4*c^4 - 42*d^4 - 72*c^2*d^2) + tan(e/2 + (f*x)/2)^4*((224*c*d^3)/3 + (8*c^3*d)/
3 + (16*c^4)/3 - (98*d^4)/3 - 40*c^2*d^2) + tan(e/2 + (f*x)/2)^2*((256*c*d^3)/3 + (16*c^3*d)/3 + (14*c^4)/3 -
(97*d^4)/3 - 44*c^2*d^2) + (80*c*d^3)/3 + (8*c^3*d)/3 + tan(e/2 + (f*x)/2)^6*(16*c*d^3 + 2*c^4 - 7*d^4 - 12*c^
2*d^2) + (4*c^4)/3 - (32*d^4)/3 + tan(e/2 + (f*x)/2)*(64*c*d^3 + 8*c^3*d + 2*c^4 - 25*d^4 - 36*c^2*d^2) - 16*c
^2*d^2)/(f*(5*a^2*tan(e/2 + (f*x)/2)^2 + 7*a^2*tan(e/2 + (f*x)/2)^3 + 7*a^2*tan(e/2 + (f*x)/2)^4 + 5*a^2*tan(e
/2 + (f*x)/2)^5 + 3*a^2*tan(e/2 + (f*x)/2)^6 + a^2*tan(e/2 + (f*x)/2)^7 + a^2 + 3*a^2*tan(e/2 + (f*x)/2)))

________________________________________________________________________________________

sympy [A]  time = 28.14, size = 8950, normalized size = 45.90 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c+d*sin(f*x+e))**4/(a+a*sin(f*x+e))**2,x)

[Out]

Piecewise((-12*c**4*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**
2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 +
f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 12*c**4*tan(e/2 + f*x/2)**5/(6*a**2*f*tan(e/2 + f*x/2)**7
 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*t
an(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 32*c**4*tan(e/2
+ f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42
*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/
2 + f*x/2) + 6*a**2*f) - 24*c**4*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2
)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2
*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 28*c**4*tan(e/2 + f*x/2)**2/(6*a**2*f*tan(e/
2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4
+ 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 12*
c**4*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/
2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**
2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 8*c**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a
**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2
+ f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 48*c**3*d*tan(e/2 + f*x/2)**5/(6*a**2*f*tan(e/2 + f*x/2
)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2
*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 16*c**3*d*ta
n(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**
5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*
tan(e/2 + f*x/2) + 6*a**2*f) - 96*c**3*d*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2
 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 +
 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 32*c**3*d*tan(e/2 + f*x/2)**2/(6*a**
2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 +
f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**
2*f) - 48*c**3*d*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*ta
n(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)
**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 16*c**3*d/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f
*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*
a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 36*c**2*d**2*f*x*tan(e/2 + f*x/2)**7/(6*
a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2
 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*
a**2*f) + 108*c**2*d**2*f*x*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6
+ 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*ta
n(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 180*c**2*d**2*f*x*tan(e/2 + f*x/2)**5/(6*a**2*f*t
an(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2
)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f)
+ 252*c**2*d**2*f*x*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**
2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 +
f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 252*c**2*d**2*f*x*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 +
 f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 4
2*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 180*c*
*2*d**2*f*x*tan(e/2 + f*x/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(
e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**
2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 108*c**2*d**2*f*x*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/2)**7
+ 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*ta
n(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 36*c**2*d**2*f*x/
(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(
e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) +
 6*a**2*f) + 72*c**2*d**2*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 +
30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(
e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 216*c**2*d**2*tan(e/2 + f*x/2)**5/(6*a**2*f*tan(e/2
 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 +
 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 240*
c**2*d**2*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/
2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2
+ 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 432*c**2*d**2*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 1
8*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e
/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 264*c**2*d**2*tan(e/
2 + f*x/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 +
42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(
e/2 + f*x/2) + 6*a**2*f) + 216*c**2*d**2*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 +
f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30
*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 96*c**2*d**2/(6*a**2*f*tan(e/2 + f*x/2)
**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*
f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 48*c*d**3*f*x
*tan(e/2 + f*x/2)**7/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2
)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2
*f*tan(e/2 + f*x/2) + 6*a**2*f) - 144*c*d**3*f*x*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f
*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x
/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 240*c*d**3*f*x*tan(e/2 + f*x
/2)**5/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2
*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f
*x/2) + 6*a**2*f) - 336*c*d**3*f*x*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x
/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a*
*2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 336*c*d**3*f*x*tan(e/2 + f*x/2)**3/(6*a**2
*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f
*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2
*f) - 240*c*d**3*f*x*tan(e/2 + f*x/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a*
*2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 +
 f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 144*c*d**3*f*x*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/
2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**
2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 48*c*d**3*f
*x/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*t
an(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2
) + 6*a**2*f) - 96*c*d**3*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 +
30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(
e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 288*c*d**3*tan(e/2 + f*x/2)**5/(6*a**2*f*tan(e/2 +
f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42
*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 448*c*d
**3*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*
x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a
**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 672*c*d**3*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*
tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/
2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 512*c*d**3*tan(e/2 + f*x/2)**
2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*ta
n(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2)
 + 6*a**2*f) - 384*c*d**3*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*
a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2
 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) - 160*c*d**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*t
an(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2
)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 21*d**4*f*x*tan(e/2 + f*x/2)**
7/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*ta
n(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2)
 + 6*a**2*f) + 63*d**4*f*x*tan(e/2 + f*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 +
 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan
(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 105*d**4*f*x*tan(e/2 + f*x/2)**5/(6*a**2*f*tan(e/2
 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 +
 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 147*
d**4*f*x*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2
 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 +
 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 147*d**4*f*x*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*
a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2
 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 105*d**4*f*x*tan(e/2 +
 f*x/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*
a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2
 + f*x/2) + 6*a**2*f) + 63*d**4*f*x*tan(e/2 + f*x/2)/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2
)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2
*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 21*d**4*f*x/(6*a**2*f*tan(e/2 + f*x/2)**7 +
18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(
e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 42*d**4*tan(e/2 + f
*x/2)**6/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a*
*2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 +
 f*x/2) + 6*a**2*f) + 126*d**4*tan(e/2 + f*x/2)**5/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)*
*6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f
*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 196*d**4*tan(e/2 + f*x/2)**4/(6*a**2*f*tan(e/2
 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 +
 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 252*
d**4*tan(e/2 + f*x/2)**3/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f
*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*
a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 194*d**4*tan(e/2 + f*x/2)**2/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*t
an(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2
)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a**2*f) + 150*d**4*tan(e/2 + f*x/2)/(6*a
**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5 + 42*a**2*f*tan(e/2
+ f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*tan(e/2 + f*x/2) + 6*a
**2*f) + 64*d**4/(6*a**2*f*tan(e/2 + f*x/2)**7 + 18*a**2*f*tan(e/2 + f*x/2)**6 + 30*a**2*f*tan(e/2 + f*x/2)**5
 + 42*a**2*f*tan(e/2 + f*x/2)**4 + 42*a**2*f*tan(e/2 + f*x/2)**3 + 30*a**2*f*tan(e/2 + f*x/2)**2 + 18*a**2*f*t
an(e/2 + f*x/2) + 6*a**2*f), Ne(f, 0)), (x*(c + d*sin(e))**4/(a*sin(e) + a)**2, True))

________________________________________________________________________________________